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Deep-inelastic structure functions: Reconstruction from Mellin moments 
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We summarize the method of calculating Mellin moments of deep-inelastic structure functions in perturbative 
QCD. We briefly discuss all steps to a complete analytical reconstruction of the perturbative corrrections in 
x-space. 



1. INTRODUCTION 

The precision measurement of structure func- 
tions and their scale evolution in deep inelastic 
scattering |l]-|j| is of basic importance for testing 
perturbative QCD. By now, experimental data al- 
lows for an accurate determination of the parton 
distribution functions, which serve as input for 
numerous other hard scattering processes. There- 
fore, in particular with respect to the upcoming 
LHC experiments, it is highly desirable to further 
reduce the theoretical uncertainties on these very 
fundamental quantities. This task requires the 
calculation of higher order perturbative QCD cor- 
rections. In order to perform a consistent next-to- 
next-to-leading order analysis, which is expected 
to significantly reduce the theoretical uncertain- 
ties Q, one has to calculate the still unknown 
anomalous dimensions of deep-inelastic structure 
functions to three loops. 

The method of calculating Mellin moments of 
deep-inelastic structure functions is directly con- 
nected with the early days of QCD G|-0|. Since 
then, it has been applied several times to deter- 
mine higher order perturbative corrections |p[-p^| , 
some of them related to sum rules. The method 
rests on the operator product expansion (OPE), 
which allows the calculation of either all even or 
all odd integer moments N of the anomalous di- 
mensions and coefficient functions. This informa- 
tion by itself then uniquely determines the quan- 
tity under consideration for general N in Mellin 
space or, equivalently its corresponding expres- 
sion in x-space. 

In the following, we will summarize the calcula- 
tion of the anomalous dimensions and coefficient 



functions in terms of harmonic sums ||16|.|l7| and 
the subsequent analytical reconstruction based 
on an inverse Mellin transformation to harmonic 
polylogarithms |18[ in a;-space. This procedure as 
detailed in ]l9| , |20| has recently been used |li| to 
check the original calculation of the two-loop co- 
efficient functions pl[-p4|, which was performed 
with conventional techniques. Currently, work is 
in progress to obtain the anomalous dimensions 
and coefficient functions at three loops by this 
method. 

The present approach of analytical reconstruc- 
tion nicely complements the use of Bernstein 
polynomials [^5j for a direct fit of Mellin moments 
to experimental data. The latter method is par- 
ticularly useful if only numerical results for a fi- 
nite number of fixed moments are available, as it 
is presently the case |H| . 

2. METHOD 

Deep-inelastic structure functions are defined 
by the commutator of two local electroweak cur- 
rents j(y) and j(z), sandwiched between hadronic 
states and Fourier transformed into momentum 



space. In the Bjorken limit, for Q 2 



and 
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x fixed, the OPE allows to express this current 
product in an asymptotic expansion around the 
lightcone (y — z) 2 ~ into a series of local com- 
posite flavour non-singlet quark O a , singlet quark 
O q and gluon operators O s of leading twist and 
spin N. The same OPE also holds for the forward 
Compton amplitude T of boson-hadron scatter- 
ing. Standard perturbation theory applies to T, 
which can be written as a series in 1/x in the 
Euclidean region. Then, by means of the optical 
theorem one obtains a direct relation between the 
parameters of the OPE and the Mellin moments 



of the structure functions. For F2 we can write 



F 2 N (Q 2 ) = J dxx N - 2 F 2 {x,Q 2 ) 



(1) 



j=">q,g 



and similar relations define and i^. 



Here 



C 2 j denote the coefficient functions and A 3 p N the 
spin averaged hadronic matrix elements of the op- 
erators O a , O q and O 6 . 

In the dispersive approach the derivation of 
eq.(Q) uses symmetry properties of the Compton 
amplitude T under exchange x — > —x. As a con- 
sequence, eq.(Q) is restricted to only either the 
even or the odd Mellin moments of F 2 . By ana- 
lytic continuation on the other hand, all moments 
in the complex TV-plane are fixed if the infinite set 
of either even or odd moments is known. To be 
precise, for unpolarized electron-proton scatter- 
ing eq.(^) determines the even moments of F 2 F , 
while for neutrino-proton scattering the odd mo- 



ments of F^'^ and F^ +vr and the even mo- 
ments of F^ p+vp and Ff^ are fixed. 

The coefficient functions and the renormalized 
operator matrix elements in eq.(Q) both satisfy 
renormalization group equations. Due to cur- 
rent conservation they are governed by the same 
anomalous dimensions, 
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Ojk 



0, 



(3) 



where j = a, q, g, while (3 and jjk represent the 
QCD /^-function and the anomalous dimensions. 
Both are calculable order by order in a s in per- 
turbative QCD as well as the coefficient functions 

2,j ■ 

The anomalous dimensions jjk in eqs.(||), (||) 
determine the scale evolution of deep-inelastic 



structure functions and, dependent on the par- 
ticular scattering process, one considers opera- 
tor matrix elements corresponding to the quark 
flavour non-singlet and singlet distributions qj s , 
q„ s and q s . Their scale evolution is governed by 
linear combinations, 



V I v 

7qq ^ 7qq 



7qq 7qq + n f (iqq 7qq) 
7qq + 7qq + «/ (7qq + 7qq) 



(4) 
(5) 
(6) 



and the flavour singlet quark distribution q s 
mixes with the gluon distribution. 

The actual calculation of the anomalous dimen- 
sions "fjk and the coefficient functions C 2 j is per- 
formed in dimensionally regularized perturbation 
theory with the external hadron states replaced 
by quark and gluon states of momentum p. By 
means of the OPE, the forward Compton ampli- 
tude for the scattering off a parton can be ex- 
pressed in terms of the coefficient functions of 
eq. (Q) and renormalized partonic operator matrix 
elements A* N . To be specific, we write for the 
flavour singlet case, 



T2 }P (x,Q 2 ,a s ,e) = 

E E (i)" c Z ($•■"•' 

N j,fc=q,g V 7 



(7) 



( a «i^) A k p N (a s ,fi 2 ,e) +0{p 2 ) 



where p = q, g and 0(p 2 ) denotes higher twist 
contributions. The Z^ represent the matrix of 
renormalization factors which contain all poles 
in 1/e. They relate bare to renormalized oper- 
ators according to O barc = Z O rcn and provide 
the anomalous dimensions via, 
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Z Z~ 



(8) 



The gauge invariant operators O q and O s in 
eq.(^|) mix under renormalization with unphysi- 
cal operators |^,^|, which are BRST variations 
of some operators or else disappear by use of the 
equations of motion. However, if one calculates 
quantities related to physical 5-matrix elements 
with physical polarization and on-shcll momenta 



such as in eq. (0) , these unphysical operators van- 
ish and have therefore been omitted in eq.(Q). 

The coefficient function C^j and the renor- 
malization factors Z^ k are calculated from eq.(|7]) 
using the method of projection [^,^9|. This 
method relies on dimensional regularization and 
amounts to the application of the projection op- 
erator, 
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N 
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d 
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where q^ 1 
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dpf* 1 ■ ■ ■ dp^>> 



(9) 



p=0 



'} is symmetrical and traceless 



part of the tensor q^ 1 • ■ ■ q^ lN . 

It is obvious, that on the right hand side 
of eq.(0) the AT-th order differentiation in Vn 
projects on the precisely the iV-th moment, which 
is the coefficient of l/(2x) N being (p-q/Q 2 ) N . All 
higher powers of p -q/Q 2 vanish after nullifying 
the momentum p. Furthermore, Vn does not act 
on the renormalization constants Z^ k and the co- 
efficient functions being only functions of N, a s 
and e. However, the nullification of the momen- 
tum p in Vn effects the partonic matrix elements 
A k N . There, it eliminates all diagrams contain- 
ing loops, as they become massless tadpole dia- 
grams, which are put to zero in dimensional reg- 
ularization. Only the tree level operator matrix 
elements A 1 ^^ survive. 

Additionally in the flavour singlet case this de- 
couples the operator mixing. Thus, eq.(|7]) pro- 
vides two independent identities to separately de- 
termine Z™, Z^ and Z<K zss. We should be 
aware however, that an Hoop calculation eq.(0) 
alone does not suffice to give the full information 
about the renormalization constants Z gq and Z gs . 
They are determined only up to (/ — l)-loops. This 
limitation is due to the gluonic coefficient func- 
tion C-fg being zero at tree level, as photons cou- 
ple directly to quarks only. To resolve this situa- 
tion and to extract the anomalous dimension 7 gq 
and 7 gg to /-loops we also calculate Green's func- 
tions in which the photon is replaced by an ex- 
ternal scalar particle <j> that couples directly only 
to gluons |15|. These Green's functions are ob- 



tained from a gauge invariant interaction term 
<t> F Zv F a V in thc Q CD Lagrangian, F« v being the 
QGD field strength. They can be expressed in 



partonic invariants T^ jP and satisfy an OPE sim- 
ilar to eqs.((l|),(0) with the same singlet operators 
O q and O g but with different coefficient functions 
C^p, where the gluonic coefficient function C^ g 
starts now already tree level. Thus the invariants 
T^p provide us with the necessary renormaliza- 
tion constants Z m and Z se of the singlet oper- 
ators to /-loops and complete the determination 
of all anomalous dimensions in the flavour singlet 
case to the desired order. 

Use of the method of projection to calculate all 
Feynman diagrams contributing to the forward 
partonic Compton amplitude, that is to say the 
application of the projection operator Vn to the 
integrands on the left hand side of eq.(^) pro- 
ceeds as follows. All ultraviolet divergences as 
well as those of infrared origin as the momen- 
tum p — > are dimensionally regularized. At the 
same time, the nullification of p drastically simpli- 
fies the topological structure of the corresponding 
graphs, which is reduced to the type of self-energy 
diagrams. Although the latter appear with sym- 
bolic powers of scalar products in the numerator 
and denominator, they are by far easier to calcu- 
late. 

In practice Jl^^^l this is performed by means of 
integration-by-parts identities and scaling identi- 
ties for the N-th Mellin moment of the diagram 
to be calculated. Combined these identities give 
rise to difference equations, which relate the N-th 
Mellin moment of the diagram under considera- 
tion to lower Mellin moments and diagrams of a 
simpler topology. Thereby they introduce a nat- 
ural hierarchy in the set of diagrams. All differ- 
ence equations can be summed recursively and 
the answer is expressible in terms of a basis of 
harmonic sums of a given weight. The 

whole approach exploits powerful summation al- 
gorithms [ fl6| for a large class of single parameter 
sums which are reducible to the basis of harmonic 
sums. 

As a demonstration of the method discussed 
thus far and as a check on l2l]-[24[|, we have 
performed a recalculation of the two-loop coef- 
ficient functions. The contributing Feynman di- 
agrams have been generated with QGRAF |30| 
and all recursion relations for the evaluation of 
the individual topologies have been programmed 



in FORM The nested sums which one en- 

counters in this way are solved with the SUM- 
MER algorithm in terms of the basis of har- 
monic sums. In addition, there is the possibility 
to perform checks at all stages of the calculation 
by means of the standard MINCER routine H) . 
This last feature is very important from a prac- 
tical point of view as the debugging is greatly 
simplified. 

To choose a specific example, we consider the 
case of initial state quarks in D — 4 — 2e dimen- 
sions and calculate with the help of eqs.(|^),(^) the 
projected invariant 7-VT2 jq . After performing the 
coupling constant renormalization one is only left 
with collinear singularities associated with the 
initial state quarks. Expanding Pjv^.q in a s and 
e gives at leading order 



At first order in a s 
e and find, 

(i) _ 



(10) 

we have to expand up to order 



T, 



2,q 



(11) 
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where the factor S e — exp(e{ln(47r) — 7e}) and 
7 q q\ are the coefficients of (a s /4ir) n in the 
expansion of the anomalous dimension 7 qq and 
the quark coefficient function C^ q . 

At second order in a s , we need to split up the 
contributions into flavour non-singlet and singlet 
parts, 



T, 



;2),s 

2,q 



= T, 



(2),ns,+ 



2,q 



rrn(2),pS 

J 2,q • 



(12) 



Allowing for electroweak interactions, one can 
consider in the non-singlet case the structure 
functions of different physical processes, i 7 ^ p±T7P , 
which implies a distinction of even and odd mo- 



ments, and if ) ' ns '" 



have 

(2),ns,± 



'2,q 



to be precise. We 
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where 7^q' V denotes the combinations 7^ ± 7^ 
as given in eqs.(Q)-(g). 

The pure-singlet contributions on the other 
hand give at second order in a s , 



T?): ps = n f 



L 2,q 



( a sY Q 2 I M 



2 c 



(14) 



+c. 



> 'qg 'gq ' 
(2),ps (o) (l) 

2,q ' 7gq «2,g 



y (o) c (i) 



and represents the combinations 7q q ±7qq of 



eqs.(|,(l). 

With eqs.©,© and at hand, we 

have at this point calculated all anomalous di- 
mension and coefficient functions up to two loops 
in Mellin space as analytical functions of the mo- 
ment N. Depending on the physical process un- 
der consideration, the quantities of interest are 
defined for either even or odd N. In a follow- 
ing step, we perform an inverse Mellin transfor- 
mation to obtain the corresponding expressions 
in x-space. As discussed in ]l8| p3|, we can do 
this step analytically as there exists a one-to-one 
mapping p^j2p[ | of harmonic sums of weight w + 1 
to the set of harmonic polylogarithms of weight 
w multiplied by 1/(1 ± x). The latter one is the 
natural class of functions for calculations of deep- 
inelastic structure functions. Since the resulting 
expressions in x-space are unique, we may subse- 
quently execute another Mellin transformation to 
go back to A-space, which then serves as an an- 
alytical continuation to all non-negative integer 
values in N . 

In particular, we are then ready to determine 
individually all anomalous dimensions which con- 
trol the scale evolution of the structure func- 
tions in the flavour non-singlet case as entering in 
eqs.(||)-(||). This is easily done now, as we obtain 
from the even moments of F^ p the sums 7l}j v and 

z ^ 1 qq 

7 q t q ' S - The odd moments of i 7 ^ p ~ lyP on the other 
hand determine the differences 7^q' V and 7qq , the 
latter being zero up to two loops. As empha- 
sized above, the analytical continuation provides 
us with expressions valid for all non-negative in- 
teger values of N, such that at this point we can 
simply take the sums and the differences of 7 ( j q ' V ' S 




and 7 qq ' V ' S to obtain all quantities T^qiTqqiTqq 
and 7 |qOfeq S .(@)-(D. 

Similar reconstructions in Mellin space and x- 
space may be carried out for the two loop coef- 
ficient functions in eqs.(|l^),(|l4|). To summarize, 
we find complete agreement with the published 
results for the anomalous dimensions 
and for the coefficient functions 



3. CONCLUSION 

The calculation of higher order pertubative cor- 
rections to deep-inelastic structure functions is 
vitally important to improve quantitative predic- 
tions for many hard scattering processes in QCD. 

In the past, this has been done either in x- 
space or in Mellin space. Using the method of 
projection, which enables us to obtain anomalous 
dimensions and coefficient functions at the same 
time, we perform the calculation in Mellin space. 
Due to new insight into the mathematical prop- 
erties of harmonic sums and their interplay with 
harmonic poly logarithms, we are first of all able 
to solve all nested sums in Mellin space. Sec- 
ondly we can reconstruct the complete analytical 
expressions of the results in a;-space by means of 
an inverse Mellin transformation. Thereby we ex- 
ploit all advantages of working in Mellin space, 
such as the use of difference equations in the 
Mellin moment N to solve the resulting loop inte- 
grals and the possibility for independent checking 
with the MINCER routine. 

Finally, we would like to remark, that the 
method presented allows for a direct application 
to the calculation of the three loop anomalous di- 
mensions, as well as a generalization to polarized 
deep inelastic scattering 35 1. 



REFERENCES 

1. D. J. Gross and F. Wilczek, Phys. Rev. D8, 
3633 (1973). 

2. D. J. Gross and F. Wilczek, Phys. Rev. D9, 
980 (1974). 

3. H. D. Politzer, Phys. Rev. Lett. 30, 1346 
(1973). 

4. W. L. van Neerven and A. Vogt, Nucl. Phys. 
B568, 263 (2000), |hep-ph/9907472] . 



5. E. G. Floratos, D. A. Ross, and C. T. Sachra- 
jda, Nucl. Phys. B129, 66 (1977). 

6. E. G. Floratos, D. A. Ross, and C. T. Sachra- 
jda, Nucl. Phys. B139, 545 (1978). 

7. E. G. Floratos, D. A. Ross, and C. T. Sachra- 
jda, Nucl. Phys. B152, 493 (1979). 

8. A. Gonzalez-Arroyo, C. Lopez, and F. J. Yn- 
durain, Nucl. Phys. B153, 161 (1979). 

9. A. Gonzalez-Arroyo and C. Lopez, Nucl. 
Phys. B166, 429 (1980). 

10. C. Lopez and F. J. Yndurain, Nucl. Phys. 
B183, 157 (1981). 

11. D. I. Kazakovand A. V. Kotikov, Nucl. Phys. 
B307, 721 (1988). 

12. D. I. Kazakov and A. V. Kotikov, Nucl. Phys. 
B345, 299 (1990), Erratum. 

13. S. A. Larin, F. V. Tkachev, and J. A. M. Ver- 
maseren, Phys. Rev. Lett. 66, 862 (1991). 

14. S. A. Larin and J. A. M. Vermaseren, Phys. 
Lett. B259, 345 (1991). 

15. S. A. Larin, P. Nogueira, T. van Ritbergen, 
and J. A. M. Vermaseren, Nucl. Phys. B492, 



16 



17 



19 



20 



25 



26 



338 (1997), |hcp-ph/9605317 



J. A. M. Ver maseren, Int. J. M od. Phys. A14, 
2037 (1999), |hep-ph/980628C . 
J. Blumlein an d S. Kurth, Phy s. Rev. D60, 
014018 (1999), |hep-ph/9810241 . 
18. E. Remiddi and J. A. M. Vermaseren, (1999), 



Phys. 



hep-ph/9905237 . 

S. Moch and J. A. M. Vermaseren, Nu cl 
Phys. B573, 853 (2000), |hcp-ph/9912"355| . 
J. A. M. Verma seren and S. Moch, (2000) 
hep-ph/0004235| , these proceedings. 

21. W. L. van Neerven and E. B. Zijlstra 
Lett. B272, 127 (1991). 

22. E. B. Zijlstra and W. L. van Neerven, Phys 
Lett. B273, 476 (1991). 

23. E. B. Zijlstra and W. L. van Neerven, Nucl 
Phys. B383, 525 (1992). 

24. E. B. Zijlstra and W. L. van Neerven, Phys 
Lett. B297, 377 (1992). 
J. Santiago and F . J. Yndurain, N ucl. Phys 
B563, 45 (1999), lhcp-ph/9904344 . 
R. Hamberg and W. L. van Neerven 
Phys. B379, 143 (1992). 

27. Y. Matiounine, J. Smith, and W. L. van Neer- 
Phys . Rev. D57, 6701 (1998), [hep" 



Nucl. 



ph/9801224 



28. S. G. Gorishnii, S. A. Larin, and F. V. 
Tkachev, Phys. Lett. 124B, 217 (1983). 

29. S. G. Gorishnii and S. A. Larin, Nucl. Phys. 
B283, 452 (1987). 

30. P. Nogueira, J. Comput. Phys. 105, 279 
(1993). 

31. J. A. M. Vermaseren, Symbolic Manipulation 
with FORM (Computer Algebra Nederland, 
Kruislaan 413, 1098 SJ Amsterdam, 1991), 
ISBN 90-74116-01-9. 

32. S. A. Larin, F. V. Tkachev, and J. A. M. Ver- 
maseren, NIKHEF-H-91-18. 

33. G. Curci, W. Furmanski, and R. Petronzio, 
Nucl. Phys. B175, 27 (1980). 

34. W. Furmanski and R. Petronzio, Phys. Lett. 
97B, 437 (1980). 

35. S. Moch and J. A. M. Vermaseren, in prepa- 
ration. 



